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ABSTRACT 

n Carinae may be the most massive and luminous star in the Galaxy and is suspected 
to be a massive, colliding wind binary system. The CHANDRA X-ray observatory has 
obtained a calibrated, high-resolution X-ray spectrum of the star uncontaminated by 
the nearby extended soft X-ray emisssion. Our 89 ksec CHANDRA observation with the 
High Energy Transmission Grating Spectrometer (HETGS) shows that the hot gas near 
the star is non-isothermal. The temperature distribution may represent the emission 
on either side of the colliding wind bow shock, effectively “resolving” the shock. If so, 
the pre-shock wind velocities are ~ 700 and ~ 1800 km s" 1 in our analysis, and these 
velocities may be interpreted as the terminal velocities of the winds from t] Cannae and 
from the hidden companion star. The forbidden-to-intercombination (//*) line ratios 
for the He-like ions of S, Si and Fe are large, indicating that the line forming region lies 
far from the stellar photosphere. The iron fluorescent line at 1.93A, first detected by 
ASCA is clearly resolved from the thermal iron line in the CHANDRA grating spectrum. 
The Fe fluorescent line is weaker in our CHANDRA observation than in any of the 
ASCA spectra. The CHANDRA observation also provides the first high-time resolution 
lightcurve of the uncontaminated stellar X-ray emission from r, Carinae and shows that 
there is no significant, coherent variability during the CHANDRA observation. The f? 
Carinae CHANDRA grating spectrum is unlike recently published X-ray grating spectra 
of single massive stars in significant ways and is generally consistent with colliding wind 
emission in a massive binary. 
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1. Introduction 

The superluminous star r, Carinae has been observed by nearly every X-ray satellite flown 
EINSTEIN Observations (Seward et al. 1979; Seward and Chlebowski 1982; Chlebowsk, et a ) 

SI the star’s X-ray emission from the rest of the Carina Nebula and mapped on pent Me 
emission centered on the star and an elliptical X-ray ring around the star extending ou ‘ 

GMGA (Koyama et al. 1990) and BBXRT provided the first clear measures of the strong Fe me 
^Hve of thermal emission probably produced by shocked gas. ROSAT discovered the van able 
nature of the hard source (Corcoran et al. 1995), while ASCA discovered > 50x solar abundance 
of nitrogen (Tsuboi et al. 1997) in the outer homunculus and fluorescent Fe Kern, ss, on (Corcoran 

et al 1998) unresolved from the strong thermal line. A 2-10 keV X-ray >6 curv « ° a ' ne 
RXTE (consisting typically of 1 observation per week since 1996) revealed small amplitude period 
^fj-th 7 = 85 ays (Corcoran et al. 1997) and confirmed that the variability seen in 992 by 
ROSAT recurred (Ishibashi et al. 1999; Corcoran et al. 2001) on the same 5.5 year period which fits 
the He 1 0830A line variability reported by Damineli (1996). These i recent X-ray - J— 
along with variations in the radio (Duncan et al. 1995) and near-IR pammd, et 200 ) all 
i th.t n Carinae mav be the Galaxy’s most massive binary system (Damineh, Conti and 
LmmT 1997' Corcoran et al.‘2001) in which the hard X-ray emission is produced by shocked gas in 
the region where the wind from , Carinae collides with the wind from the 

observation of n Carinae by the Advanced CCD Imaging Spectrometer (ACIS) on CHANDRA n 

September 1999 (Seward et al. 2001) provided the first image of the X-ray regions a 

JT T” though' due to degradation produced by charged partide damage of the CCDs, tins 

observation could not be precisely calibrated spectrally. 

Here we report the first calibrated observation of r) Carinae with the 
Grating Spectrometer (HETGS, Canizares et al. 2001, in preparation on the CHANDRA X . ray 
observatory (Weisskopf, O’Dell, k van Speybroeck 1996) using the ACIS linear array (ACIS-S to 
read out the dispersed spectrum. The CHANDRA HETGS is well suited for observing the point-like 
hafd emission La V Carinae for two reasons: 1) the spectral energy distribution of the emission 
L nic el y in the HETGS bandpass and 2) this emission is thought to be dominated by therma 
processes which should produce strong line emission in the dispersed spectrum. ur eep ( . 

hour) observation of rj Carinae with the HETGS provides the first spectrally-resolved measure o 
the X-ray spectrum, allowing for the first time a detailed definition of the temperature distribution 
density Id chemical abundance of the hot, unresolved emission from the star. In this firs repor 
we discuss the overall spectral morphology of the unresolved X-ray emission, examine e emission 
ta Z^,a«u,e uud density diagnostics, and provide an uninterrupted lighten, ve of the unresolved 

X-ray source over the length of the observation. 
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2. The Observation 

The CHANDRA HETGS+ACIS-S observation of r/ Carinae was performed on 19 November 
2000 - 20 November 2000. The total exposure time of the observation was 89,546 seconds The 
spacecraft roll was chosen so as to avoid contamination of the dispersed spectra y ot er rig . 
ray sources in the field. The data were cleaned and processed using the standard pipeline processing 
available at the CHANDRA X-ray Center. Images and spectra were extract rom 
events using the Chandra Interactive Analysis of Observations (CIAO) analysis package. 

The zero^-order image, color-coded by X-ray energy is shown in figure 1. This image is similar 
to that previously published (Seward et al. 2001) though the energy calibration is better in the new 
image. In particular, the ACIS-S zero‘ k -order image shows the soft elhptic^ “shell ^of emission 
surrounding the hard X-ray core, which is unresolved to ACIS at scales of ~ 0.5 ~ cm ~ 

AU. The X-ray flux in the elliptical “shell” is very non-uniform, bright in the south and west (near 
the “S-ridge” and the “W-condensations” in the outer debris field, Walborn, Blanco & ac cry 
1978) and faint in the north and east. The hard, unresolved “core” appears to be surrounded by a 
halo of X-ray emission at moderate X-ray energies. 

3. The X-ray Grating Spectrum of r) Carinae 

The ACIS-S image of the dispersed spectrum shows that the hard source is unresolved to 
CHANDRA ; in particular, there is no observed dispersed spectrum from any other source of emission 
in the ACIS-S field aside from the unresolved, hard core source associated with r, Cannae. JVe 
extracted the dispersed medium energy grating (MEG) and high energy grating (HEG) J a 

from the X-ray event file using CIAO. Both plus and minus orders for the first, second an 
order MEG and HEG source spectra were extracted, along with appropriate background spectra. 
Figure 2 shows the MEG +1 order spectrum, while figure 3 shows the HEG +1 order ,n the vicinity 
of the Fe K line. The MEG and HEG spectra of t? Carinae represent the first X-ray spectra o is 
star uncontaminated by extended soft emission. The MEG spectrum shows very little emission 
from the unresolved source at energies less than 1.5 keV (A > 8.5 A) due to strong absorption, 
reveals a significant X-ray continuum in the 1 - 8A range and the presence o strong me emission 
from lines of S XV-XVI, Si XIII-XIV, Mg XII, CaXIX, and Fe K. Strong forbidden lines of S XV, S. 
XIII and Mg XII were detected. In the HEG spectrum, the Fe line region shows a thermal emission 
line produced by Fe XXV and a blend of fluorescent Fe I Kai+Ka 2 hnes centere at .. 


3.1. Modelling the Spectrum 

We attempted to lit the MEG spectrum with a combination o( optically-th'm thermal emission 
models using the XSPEC analysis package (Arnaud 1996). Because the dispersal spectrum has 
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many bins with few counts and because background does not contribute significantly n the ene^y 
range of interest (E > 1.5 keV or A < 8 A), we used a modified version of the C-statistic (C 
1979) on the total (non-background-subtracted) spectrum instead of the X statistic. Previous 
analyses indicated that the emission at energies above 1.5 keV could be fit by emission at a single 
temperature (for eg., Corcoran et al. 2000). As shown in Figures 4 and 5, a variable abundance 
single temperature thermal model (VMEKAL-XSPEC, Mewe, Kaastra k Liedahl 1995) provide 
a gLd fit to the continuum emission and most of the resolved emission lines. The parameters of 
the best fit single component model are given in Table 1. The single component temperature ,s 
kT = 4 4 keV with a column density of N H = 4.9 X 10 22 cm' 2 , in reasonable agreement with earlier 
results However, we found that no isothermal model could simultaneously match the strengt is 
lAt Mke and BA lines of Si in .he 4. 5 A < A < 6 A range. Fitting both the He-„ e and 
H-like ions require the addition of at least one additional thermal component. On, best fit 2- 
temperatur. model is given in table 2 and shown in figures 4 and 5. Tins 2-‘eml»rat«re ». 
adequately matches the strength of both the Holike and H-hke I, no. This ,s the first tin. that a 
unresolved X-ray emission from , Carinae has been shown to require a nonosothermal tempe atur 
distribution. The maximum temperature we derive, kT as 8.7 keV, » nearly a factor of 2 large 
than most other published temperatures for q Carinae, and to on, knowledge represents the highest 
temperature ever associated with an early-type star. 

In on, modelling of the X-ray spectrum we held most elemental abundances at their sola, values 
(Anders ft Grevesse 1989), but allowed the abundance of Si, S and Fe (all of which have strong 
lines in the X-ray spectrum) to vary. The derived abundant for the 1 and 1 — ure mo^s 
fo, these 3 elements are given in table 1. In each ease the abundances of the S S, Fe e 
slightly non-solar: Si needs was found to be slightly overabundant, S significantly overabundant 
(by about 70%), while Fe was slightly underabundant. 


3.2. The He-like Lines and the f/i Ratio 

The ratio of the intensity of the forbidden component to the intercombination component of 
the He-like lines (the //• ratio) is a density diagnostic (Gabriel ft Jordan 1969), t oug e ra m 
may be also increase b, UV photoexcitation (Kahn et al. 2001) which car .suppress the 
fine and enhance the intercombination line. We fit the Si XIII, S XV and Fe XXV fines ,n XSPEC 
by first isolating the wavelength region around the line complex of interest and ustng <» 

model the lines with inclusion of a power law component to drscribe the background The fits are 
shown in figures 6-8 and the fit parameters are given in Table 2. The intercom!,, nation hues for 

all three ions are weak and our analysis only yields upper limits for the intensities of these lines. 

, / . i f . 1 o f or c; XIII > 2 0 for S XV and > 2.5 for Fe XXV. These values 
We measure f/i values of > 1.6 tor bi Alii, > * o c c vv ^ . t- \\y 

correspond to electron densities of < 10'* cm-" fo, Si XIII and < 10'= cm for S XV and Fe XXV. 
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3.3. The Fe Fluorescent Line 

We fit the Fe I fluorescent line from the HEG spectrum in XSPEC by isolating the Fe K region 
and using a single gaussian line plus power-law background component. The line parameters are 
given in Table 1. The measured equivalent width of the line is only 39 eV, which is about half the 
value of the smallest equivalent width seen by ASCA (Corcoran et al. 2000). 


4. The X-ray Lightcurve 

From the zeroth-order image on the S3 chip we extracted a lightcurve in a circular 5" diam- 
eter region centered on the unresolved source. We also extracted a background lightcurve from a 
source free circular region of diameter 17" centered just beyond the outer elliptical emission region 
surrounding the hard unresolved source. We used bins of 1000 seconds for both the source and 
background lightcurve extraction; this yielded about 200 gross source counts and about 10-20 back- 
ground counts in each bin. The background subtracted X-ray lightcurve of the unresolved source 
for the duration of the CHANDRA observation is shown in figure 9. This is the first uninterrupted 
lightcurve of the hard unresolved source spatially resolved from nearby contaminating emission. 
While there is some evidence for sporadic variability there is no evidence for coherent variations in 
the hard emission from the unresolved source during this observation. 


5* Discussion 

Recent evidence suggests that 77 Carinae may be a massive binary star, in which variable X-ray 
emission is produced by shocked gas at the interface where the wind from 77 Carinae collides with 
the wind of its companion (presumably some less massive early-type star). The line complexes 
resolved in the HETGS spectra provide unique information about the physical condition of the 
X-ray emitting region and thus on the single or binary nature of 77 Carinae. In the colliding wind 
model the wind from 77 Carinae forms a bow shock around the companion, since the companion’s 
wind is probably weaker than 77 Carinae’s wind, i.e. 77 = (MV r 0 o) c /(^V 00 )^ < 1; here M represents 
the wind mass loss rate and V x the wind terminal velocity, and the subscripts 77 and c refer to 77 
Carinae and the companion, respectively. In published models (Damineh et al. 2000; Corcoran et 
al. 2001) the mass loss rate from r, Carinae is thought to be at least a factor of 10 larger than that 
of the companion, while the terminal velocities are probably within a factor of 3 (see below), so 

that 77 < 0.3. 

The mere presence of numerous strong emission lines confirms that the X-ray emission is 
dominated by thermal processes from shocked gas. The dominance of thermal processes in the 
CHANDRA energy range suggests that there is little or no contribution from non-thermal processes 
(which could be produced via inverse compton scattering of photosphenc UV photons if Fermi 
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acceleration of the electron population in the shock is important)- The °[ 

the X-ray emission as measured by CHANDRA provides a measure of the pre-shock wind velo 

for both stars. In a colliding wind binary, the temperature of the shocked gas ^^tLle and 
kT (keV) * 2.59^V 2 (see the eq. 51 in Usov (1992)) where /x is the mean mass per particle and 

V is the stellar wind speed in 1000 km For simplicity we assume „ « U though this value 

may not be truly representative of the chemical composition of the wind of either r, Cannae or its 
companion We neglect the motion of the shock stagnation point due to the orbital motion of he 
stars, since orbital motion has a significant influence on the derived ^o^when 

stars are near periastron (and the stellar velocities are greatest), wh.le the HETGS observation was 

derived kT » 1.1 and « 8.7 keV, suggest velocities of ~ 700 km s and ~ ’ 

!f the wind from , Carinae is thought to be ~ 500 km s" (Hil.ier et al. 
with the lower velocity we derive from our X-ray analysis. A velocity of ~ 500 km s is su 
to produce the kT ~ 1 keV emission, but cannot produce the hotter emiss.on revealed in the MEG 
p^rum Naively, we interpret the low-temperature component as shocked emission from he 
relatively slow wind of , Carina., and the high-temperature component as shoe ed «— by the 
faster wind from the companion star, so that * 700 km s"* and Foe,, ~ 1800 km . 

The f/i ratios derived from our analysis of the He-like ions are large and show that the 
line forming region has a density < 10 14 cm" 3 and that enhancement of the intercomb, nation 
component by UV photoexcitation is unimportant. This suggests that the line forming region 
CafTflr from fife sfellar photosphere, where wied duties and UV » - 
broadly consistent with published descriptions of the orb.t (Dammeh et al. 2000, Corcoran . 
2001) 'in which the shocked emission forms fa, from either star (though closer to the companion 
Z J “ Carinae) in a region of low particle and UV photon density. In the colhdmg w.nd model, 
1 apex of the bow shock is located a. a distance of d/(l + y/h) f™" V C “ ina * d * ^ 
L a 7a«on of the two stars. At the phase of the CHANDRA observation, * = ™ 

are separated by about 4 X 10- cm, using the orbital elements of Conroran e, ^ (2001)^ At th.s 
time the bow shock is about 2.5 X 10 14 cm from V Cannae, and about 1.5 X 
companion. The density of the wind at a distance r from the photosphere is » = _M/(4 7 rr V M m„l 
where M is the mass loss rate and V the wind velocity. Using M n - 10 M @ y °o,n 

km s'* (Hillier et al. 2001), then the density of the wind from „ Cannae at ^ of . _ 
cm is only n « 8 X 10 8 cm" 3 . Using M c = 10 “ 5 M 0 yr" 1 (Corcoran et al. 2001) and W - 1800 
km s- 1 as appropriate values for the companion’s wind, the density of the companion s win a a 
distsnee r =1 5 X 10- cm is also only about 10» cm' 3 . The predicted dens, ties are cons, stent with 
the upper limits derived from ou, analysis of the f/i line intensity rat, os. However, ,f the actua 
density of the emission region is near the limit implied by the f/i ratios, th.s might indicate 
the assumed wind momentum balance is not quite right. 

In contrast, newly published X-ray grating spectra of 6 1 Ori C (Schulz et al. 2000), ( On 
(Waldron & Cassinelli 2001) and < Pup (Kahn et al. 2001) all have f/i ratios which are lower 
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the values we derive from the i, Carinae X-ray line spectrum. None of these stars are known to show 
any colliding wind effects, and apparently, in these massive stars, the X-ray lines form relative y 
near the stellar photosphere, probably within a few stellar radii. This is not the case for rj Carinae. 

The strength of the Fe fluorescent line is related to the column density of scattering material 
bv EW » 2 3iV 24 keV, where EW is the equivalent width of the line in keV and N 2 a the total 
column density of cold material in units of 10* cm- (Kal.man 1995). The ASCA spectra suggest 
that EW « 4-7A r 24 for r/ Carinae outside of eclipse (Corcoran et al. 2000). Thus at the time of the 

HETGS observation the equivalent width of the Fe fluorescent line implies N 24 « 0.01 - 0.03, i.e 

. r i j ic IVtt '“-'I 3 x 10 22 cm 2 . This value is in fair agreement 

that the column density of cold material is Ah ^ 

with the value of N„ = 5 X 10 22 derived from fitting the X-ray continuum in the MEG spectrum, 
suggesting that some of the same material which produces the X-ray absorption also is responsi e 
for producing the Fe fluorescent emission. The equivalent width of the line which we derive here 
is much smaller than any values given in either Corcoran et al. (2000) or Corcoran et al (1998). 
This may be the result of the difficulty in determining the actual width of the fluorescen me in 
the ASCA spectra since the line is not resolved from the thermal component. However if this 
variation is real, it may represent a real decrease in the scattering optical depth. Such a decrease 
is not unexpected in the binary model since at the time of the HETGS observation the companion 
is nearly in front. Since the companion’s wind is less dense this could produce a decrease in the 
scattering optical depth at the time of the HETGS observation. 


6. Conclusions 

We have presented here the first high-energy X-ray grating spectrum of the supermassive star 
r, Carinae. This grating spectrum shows strong line emission from H-hke and He-hke i ions of S 
Si Mg, Ca, and Fe, and confirms the presence of the Fe fluorescent line discovered by - ASCA. The 
forbidden lines of the He-like ions are strong and the intercomb.nation lines weak, indicating tha 
the line forming region is far from the stellar photosphere. These results are all consistent with the 
current picture of r, Carinae as a colliding wind binary. We expect that interesting variations in 
the X-ray spectrum will occur as the two stars approach periastron (which should next occur on 
June 20, 2003). In particular, we expect that the strength of the forbidden lines will weaken near 
periastron, as the density in the X-ray region increases and as the UV photosphenc flux in the line 
region intensifies. Additional HETGS observations as the stars approach periastron, especially in 
conjunction with UV spectroscopy, will provide unique and significant constraints on the orbital 
geometry, the strengths of the winds from both stars, and the evolutionary stages of the stellar 

components. 


We would like to thank the CHANDRA staff at the SAO/CfA and the HETGS team at MIT 

for their help in scheduling the observation and in analyzing the data. We made use of the F O s 
suite of software supported by the HEASARC, and of the ASCA data archives also maintained at 
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Fig 1.- Zeroth order image, color coded for X-ray energy (red - low energy; green -moderate 
energy; bine - high energy). The diagonal band at the top of the image rs the dispersed med.um 

energy grating (MEG) spectrum. 

Fig . 2 ._ +1 order Medium Energy Grating (MEG) apcctrum of „ Carinae. Strong lines from H- 
and H e-like ions of Mg, Si, S, Ca and Fe are indicated. 

Fie 3 - Iron K line region from the +1 order High Energy Grating (HEG) spectrum of f Carinae. 
Emission is seen from He-like iron originating in gas a, a temperature of ~ 8 keV. In add,,, on 
fluorescent line emission from cool iron is also clearly detected. 

Fie 4 - Fits to the MEG +1 order spectrum. A two-temperature fit using a variable abundance 
cl onal^y o iced plasma model (VMEKAL-XSPEC, Mewe, Kaastra Sr Liedahl 1995 » shown ,» 
brnf and a single-temperature fit with the same emission model is in red. The Single temperature 
mluanno; rtodu« simultaneously the strength of the S XV, Si XII, and Si XIV hues ,n the 

range 4.5-7 .0A. 


Fig. 5 — Same as in Figure 4a, but emphasizing the S & Si line region. 

Fig. 6.— Fit to the /, t, r lines of Si XIII. 

Fig. 7.— Fit to the /, i, r lines of S XV. 

Fig. 8. — Fit to the /, i, r lines of Fe XXV (near 1.85A) and the Fe I fluorescent line (at 1.93A). 

Fig . 9 _ The background-corrected lightcurve of the unresolved hard source. 
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Parameter ________ 

~Nh 

kT 

Emission Measure 1 
Lx (2-10, absorbed) 1 
Lx (2-10, unabsorbed) 1 
Si/Si© 

S/S© 

Fe/Fe© 

m 

Emission Measure 1 2 

C —statistic 

1 calculated assuming a 
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Units 

1 Component Fit 

2 Component 

10 22 cm" 2 

4.9 

5.1(±0.4) 

keV 

4.4 

1.1(±0.1) 

10 57 cm -3 

4.0 

1.8(±0.1) 

ergs s _1 

2.6 

3.0 

10 34 ergs s _1 

3.9 

4.5 

1.5 

1.1(±0.2) 


1.4 

1.7(±0.4) 


0.6 

0.9(±0.2) 

keV 


8.7(±1.7) 

10 57 cm -3 


2.9 (±0.3) 


2035 

1619 


distance of 2100 pc (Corcoran et al. 2001) 
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Line E (keV) 


SiXIII (r) 

1.864 

SiXIII (i) 

1.853 

SiXIII (f) 

1.839 

SXV (r) 

2.460 

SXV (i) 

2.450 

SXV (f) 

2.431 

FeXXV (r) 

6.700 

FeXXV (i) 

6.682 

FeXXV (f) 

6.636 

Fel 

6.424 


Measured E (keV) 


1.866 

1.859 

1.839 

2.460 

2.450 

2.434 

6.700 

6.680 

6.640 

6.407 


Eq. Width (eV) 


38 

< 4.3 

21.4 
33.1 

< 6.0 
25.6 
261 
< 42 

42 

39.4 


Intensity 
(photons s" 1 cm" 

0.017“ 

< 0.003 
0.004 
0.021 

< 0.003 
0.006 
0.230 

< 0.050 
0.100 
0.055 


keV" 1 ) 



